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ABSTRACT 

Type Iln and related supernovae show evidence for an interaction with a 
dense circumstellar medium that produces most of the supernova luminosity. X- 
ray emission from shock heated gas is crucial for the energetics of the interaction 
and can provide diagnostics on the shock interaction. Provided that the shock 
is at an optical depth r w < c/v s in the wind, where c is the speed of light 
and v s is the shock velocity, a viscous shock is expected that heats the gas 
to a high temperature. For r w > 1, the shock wave is in the cooling regime; 
inverse Compton cooling dominates bremsstrahlung at higher densities and shock 
velocities. Although r w > 1, the optical depth through the emission zone is < 1 
so that inverse Compton effects do not give rise to significant X-ray emission. The 
electrons may not reach energy equipartition with the protons at higher shock 
velocities. As X-rays move out through the cool wind, the higher energy photons 
are lost to Compton degradation. If bremsstrahlung dominates the cooling and 
Compton losses are small, the energetic radiation can completely photoionize the 
preshock gas. However, inverse Compton cooling in the hot region and Compton 
degradation in the wind reduce the ionizing flux, so that complete photoionization 
is not obtained and photoabsorption by the wind further reduces the escaping 
X-ray flux. We conjecture that the combination of these effects led to the low 
observed X-ray flux from the optically luminous SN 2006gy. 

Subject headings: circumstellar matter — shock waves — supernovae: general — 
supernovae: individual (SN 2006gy) — X-rays: general 



1. INTRODUCTION 



There is increasing evidence for supernova shock waves propagating in dense, optically 
thick, mass loss regions. Type Iln supernovae, which have narrow lines of hydrogen and 
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other species in their spectra (ISchlegellll990f). typic ally have light curves where circumstellar 
interaction provides the power (e.g., IChugail Il992[ ). The lines are indicative of continuing 
circumstellar interaction and the luminosity implies a high circumstellar density. Narrow H 
lines with broad, symmetric wings observed at early times can be interpreted a s the result 
of electron scattering in a medium with optical depth of 3 — 4 ( IChugail 120011 ) . Type Iln 
spectral features have been observed in num ber of highly luminous supernovae, including SN 
2006gy JSmith et al.ll2007t bfek et al.ll2007f ). SN 2006gy, which was especially luminous in 
the optical, had an X-ray luminosity (or upper limit) th at was orders of magnitude below its 



optical luminosity fjSmith et al.l 120071 ; lOfek et al.l 120071 ); the question arises of whether the 



lack o f a high X-ray luminosity is consistent with strong circumstellar interaction (IKatz et al. 
201lh . 



There is the potential for X-ray emission provided that there is a viscous shock wave 
that heats the gas in the circumstellar interaction. A viscous shock front is expected to 
form in a dense wind pr ovided the optical depth in the external medium r < c/t>, s , w here 
v s is the shock velocity Jofek et al.lhoid Ichevalier & Irwinlboill iNakar k Sarilboioh . At 
larger optical depths the shock wave is mediated by radiation pr ovided that th e ratio of 
radiation to matter pressure is > 4.45 in the downstream region (jWeaverl Il976l ). which is 
the case for the shocks considered here. There is thus the pote ntial for hot gas and its 
emission from shock waves in moderately optically thick regions (IKatz et al.l 1201 ll ). Here 
we consider the X-ray emission from shocks in dense circumstellar regions, concentrating on 
the optically thick case. Any emitted X-ray emission has the possibility of being scattered 
or absorbed, changing both the supernova surroundings and the escaping X-ray radiation. 
An understanding of how the shock power is eventually radiated is crucial for understanding 
multiwavelength observations of these events. 

We discuss the shock structure and X-ray emission in Section 2 and the implications for 
observations in Section 3. 



2. SHOCK PROPERTIES AND X-RAY EMISSION 

We assume that the mass loss can be described by a steady wind flow; the actual mass 
loss is unlikely to be steady, but this case should illustrate the basic physical situation. If 
the mass loss at rate M is in a steady wind at a velocity v w , the density p w = M/4irr 2 v w = 
Dr~ 2 can be specified by a density parameter, D*, scaled to a M = 10~ 2 M yr _1 and 
v w = 10 km s -1 wind so that p w = 5.0 x 10 16 D*r~ 2 in cgs units. The optical depth in the 
wind to the dense shell position R is r w — 1.7 x 10 16 kD*/R where k is the opacity k in units 
of 0.34 cm 2 g _1 , appropriate for electron scattering with riH e / n H = 0.1. If the expansion of R 
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can b e expressed R ~ t m where t is the age (I Chevalier fc Franssonl 120031 ; IChevalier fc Irwin 



201ll ). we have v s = mR/t so that 

D* = OMVs^k^it/lO day), (1) 

where v S 4 is the expansion velocity of the shell v s in units of 10 4 km s _1 and m = 0.8 has been 
assumed. This relation with r w — 1 is shown in Fig. 1, which also shows the line r w = c/v s 
above which a viscous shock does not form. In the figure, and in the following, we assume 
k = 1. 



2.1. Radiation from the Shocked Region 

To give an est i mate of the shock velocity, we note that for the supernova model used by 
Chevalier fc Irwin! (120111 ) (p oc r -7 outer density profile), the velocity of gas at the reverse 



shock wave is 6.3 x 10 3 ££i 2 M- - 2 Z>-°- 2 (t/10 day)" - 2 km s -1 , where E 51 is the energy in units 
of 10 51 ergs and M el is the ejecta mass in units of 10 M . The velocity v of the interaction 
shell is also the shock velocity v s if the velocity of matter immediately ahead of the shock is 
negligible. There are two reasons that v s might be less than v. One is that the presupernova 
wind has a significant velocity. Wind velocities deduced from narrow Ha lin e features in 



Type Iln supernovae are generally in the range 100 — 1000 km s 1 (Table 5 of iKiewe et al. 



20121 ). Second, there may be radiative acceleration of the mass loss gas. 



To estimate the acceleration, we assume that the forward shock velocity is ~ v and that 
the shock wave cools rapidly (compared to the age). These assumptions should give the 
maximum luminosity and thus the maximal effect of radiative acceleration. The luminosity 
in this case is L c = 2irR 2 p wQ v 3 , where p W Q is the wind density immediately ahead of the 
shock wave. The acceleration is a ra( i = Lk/47[R 2 c, where c is the speed of light, so that the 
gas is accelerated to a velocity 

^ p w0 v 2 mRK 

Vrad Qradt ~ „ • (^) 

2c 



If we normalize to the shock velocity, 

V rad mDvK 771 V 

~V ~ 2Rc ~ 2"c 



t w . (3) 



The result is that v ra d ~ v if t w ~ c/v. The formation of a viscous shock requires that 
t w < c/v s , so we expect significant acceleration near the time of shock breakout, if it occurs 
in the wind region, and declining acceleration at lower optical depths. The breakout radiation 
from a hot shell of radiation can also accelerate the wind gas, but the luminosity of this event 
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i s comparable to the lum in osity that might b e generated by the continuing shock interaction 
( IChevalier &: Irwin! 120 111 ). iKatz et al.l ( 1201 ll ) obtain a similar result for the acceleration by 
breakout radiation and note that the shock is collisionless. Since the acceleration is important 
only at fairly large optical depth and we will find that X-rays do not escape at those depths, 
we neglect the acceleration here. 

The cooling processes for the postshock hot gas are expected to be bremsstrahlung 
and inverse Compton cooling. At the high densities of interest here, the dominant cooling 
emission is from the forward shock region. X -ray emission from supernova e is often discussed 
in terms the softer reverse shock emission (IChevalier &: Franssonl 120031 ). However, as the 
density increases, the reverse shock emission becomes limited by rapid cooling at a lower 
circumstellar density than at the forward shock. When both the forward and reverse shocks 
are in the radiative regime, the ratio of forward to reverse shock power is 2 (re — 3) 2 / (re, — 4), 
where re is the supernova density power law index (IChevalier fc Franssonl 120031 ). For re = 7, 
the ratio is 11. In addition, the dense shell that forms between the shock fronts can absorb 
X-ray emission produced at the reverse sho ck. In the noncooling regime , the bremsstrahlung 
emission from the forward shock region is (IChevalier fc Franssonl 120031 ) 



L b = 3 x 10 45 L> 2 (t/10 day)" 1 erg s 
In the cooling regime, the luminosity is 

L c = 3.1 x 10 44 £>^ s 3 4 erg s -1 , 



-i 



(4) 



(5) 



where v s 4 is the shock velocity in units of 10 4 km s 1 . The actual luminosity cannot be 
greater than that in the cooling case. The cooling condition is that L b > L c , or 



A, > 0.1^ 3 4 (t/10 day). 



(6) 



The cooling condition is shown in Fig. 1; it is close to where the medium is optically thick. 

Since the two sources of luminosity for inverse Compton cooling (shock breakout emis- 
sion and emission from continuing interaction) are comparable at early times and the contin- 
uing interaction eventually can dominate, here we consider only the interaction luminosity 
source. If the forward shock is cooling, the postshock luminosity is L c = 2irR 2 poVg. We as- 
sume that a fraction / of this luminosity is reradiated as approximately blackbody thermal 
radiation and take / = 0.5 as a reference value; backscattering of the radiation can cause 
deviations from this value. In the optically thick regime (t w > 1), the energy density in 
radiation is thus 

fL 



^rad 



47rr 2 c 



(7) 
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The inverse Compton energy loss per unit volume is 

k B T e 



A c = Au rad cn e a T - 



m e c z 



(8) 



where n e is the electron density in the hot gas, <jt is the Thomson cross section, ks is 
Boltzmann's constant, m e is the electron mass, and T e is the electron temperature. In the 
radiative, optically thick regime, the ratio of the bremsstrahlung cooling rate to the Compton 
cooling rate is 



^ 0.7 — 
Ac V0.5 



^sA 1~w 7 



(9) 



where equation (3-56) of ISpitzerl (119781 ) was used for A br . The boundary between the cooling 
mechanisms is shown in Fig. 1. 

Moving to the non-cooling case, we note that L c oc D in the cooling case, but the 
luminosity of the shell L = VAf, r oc D 2 in the noncooling case, where V is the emitting 
volume of the hot gas. We again assume that the supernova luminosity is primarily from 
circumstellar interaction so that the photospheric luminosity is fVA^. More specifically, we 
find 



Abr 

Ac" 



4/r 



kT 



m P c z 



2.1 ( -L 

0.5 



(10) 



The electron scattering optical depth through the hot shell, To, is the same as that through 
the preshock wind gas, t w , provided H and He are ionized and dominate the abundances, 
and the shell gas is noncooling. The value of To drops below unity in the noncooling regime. 
Inverse Compton cooling is important at higher shock velocities and wind densities (Fig. 1). 
This contrasts with the case where the radiation field is due to the supernova, independent 
of the circumstellar int eraction, when bremsstrahlung emission becomes a more important 
coolant at high density (jFranssonlll982l ) . At the higher densities in our case, inverse Compton 
again becomes important due to the strong radiation field created by the interaction. As 
indicated in Fig. 1, the condition that the shock front be cooling merges with the transition 
between inverse Compton and bremsstrahlung cooling at higher shock velocities. The reason 
for the absence of a region in which the shock is non-cooling with cooling dominated by 
inverse Compton is our assumption that the supernova luminosity is produced by the shock 
interaction. 

If the nuclei and electrons do not rapidly achieve equilibrium in the shock transition, 
the nuclei are heated to T p = 2.9 x 10 9 t>g 4 K in the shock front. The balance of electron 
cooling by inverse Compton losses with heating by Coulomb collisions leads to a temperature 
T e « 7.1 x 10 8 e^ 2 ^ 5 K, where e 7 is the fracti on of the postsho ck energy density in radiation 



and a Coulomb logarithm of 30 is assumed (IKatz et al.ll201ll ); T e may be higher if there is 
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collisionless heating of electrons. From equation (J7]), we have e 7 ~ (t w v s /c), so that when 
the viscous shock first forms e 7 ~ 1, but it declines thereafter. The electrons are heated to 
their equilibrium value (1.4 x 10 9 t>^ 4 K) when r w « 5.7iQ . Nonequilibrium is important 
only for high shock velocities (Fig. 1). The criterion for equilibrium when the cooling is 
domin ated by bremsstrahlung at lower densities is the same as that discussed by iFransson 
(119821 . equation 12) and is shown in Fig. 1. 

For t w < 1, bremsstrahlung is expected at close to the postshock temperature and, since 
the emission goes in toward the dense photosphere as well as outward, with a luminosity 
approximately equal to that from the photosphere. As the optical depth to the emission 
region increases, there is more of a chance of an outward going photon being scattered in 
toward the photosphere, so the hard X-ray luminosity declines relative to the photospheric 
emission. In the cooling case, the electron scattering optical depth through the hot, shocked 
region is r = n e <TTd, where d = {y s /4)t coo i and t coo \ = (3/2)nfcsT/A. If inverse Compton 
dominates the cooling and t w > 1, A is given by equations (J7J) and dHJ, and we have 

m P ( c \ 2 0.6 



T ° Af/j,r w m p \v s ) fr w v^ ^ 

where \i is the mean particle weight divided by the proton mass m p . We thus find that r < 1, 
so that a typical outgoing photon from the dense shell scatters at most once in the hot gas; 
since the electron energies are < m e c 2 , the photon energy is changed by a factor < 2. For 
t w > 1, photons can be scattered back through the hot region, but production of photons up 
to X-ray energies is not expected in this way because the ingoing photons are absorbed by the 
dense shell and ejecta. We thus expect that, in the regime where inverse Compton cooling 
dominates, the X-ray emission declines relative to emission at lower photon energies. This 
is also true when bremsstrahlung cooling dominates and t w > 1 because initially outgoing 
photons can be scattered back across the emission region and absorbed in the dense shell. 



2.2. Photon Interaction with the Preshock Wind 



The X-ray radiation from the hot shell must escape through the cooler unshocked enve- 
lope of circumstellar gas. Both scattering and absorption of the photons can be important. 
Comptonization affects the escape of photons at the high energy end. The maximum energy 
is (-max ~ m eC 2 /Tg S = 511 / t^ s keV, where the electron scattering optical depth r es = r w for 
our assumptions. This result is due to the fact that electron recoil gives a change of photon 
wavelength ~ h/m e c for each scattering, where h is Planck's constant, and the number of 
scatterings is ~ t^ s . Detailed calculations show tha t there is not a sharp cu t-off at the energy, 



but the spectrum becomes steep (e.g., Fig. 15 of iKylafis &: Lamblll982l ). In Section 1, we 
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noted that a viscous shock could form at r w & c/v s = SOvJ^ so that e max can be as small as 
0.57f^ 4 keV. The emission from the hot gas goes into heating and ionization of the preshock 
gas until the shock wave reaches moderate optical depth. 

Photoionization of the preshock medium is important for the absorption of X-rays from 
the shocked region. If the preshock medium is completely ionized, we expect relatively little 
absorption, while incomplete ionization leads to substantial absorption for our parameters. 
In considering the photoionization of the preshock gas, we will be assuming that the pho- 
toionization is determined by the current X-ray luminosity, i.e. that steady state conditions 
apply. This requires that the recombination time be less than the age, which is generally 
true for the high densities considered here. 

If the medium is optically thin, photoionization is related to the ionization parameter 
£ = L/nr 2 in cgs units, where n is the density and r is the distance from the luminosity source. 
For an r~ 2 density distribution, the £ parameter is independent of radius. If we are in the 
cooling regime and the ionizing luminosity is (3L C , the ionization parameter is £ ~ 1 x 10 A (3v 3 4 , 
independent of D. For lower values of D in the non-cooling regime, we have £ oc Dt~ l . The 
ionization parameter is highest in the high density regime of interest here and drops at lower 
densities in the noncooling regime. Photoionization calculations have been previously carried 
out and we briefly summarize results for an emitting gas temperature of 10 keV, or ~ 10 



K. For a value £ ~ 10 , the medi um is completely ionized ( iTarter et al.lll969t lHatchett et al. 



19761 ; iKallman &: McCrayi 119821 ) . The elements C, N, and O are complet ely ionized for 
£ > 100; ionization of the heavier elements (S and Fe) requires £ > 10 3 (lHatchett et al. 



19761). 



To ext end these results to a hotter radiation field, we used the photoionization code 



CLOUDY (IFerland et al.lll998T ). We started by using a radiation field set by a T = 10 8 K 
bremsstrahlung spectrum; the results were consistent with those described above. We went 
to higher bremsstrahlung temperatures to allow for higher shock velocities; at v s 4 = 1, the 
postshock temperature is ~ 10 9 K. At £ = 1000, changing the bremsstrahlung temperature 
does not have much effect on ionization structure or T e . For £ < 1000, a higher temperature 
luminosity leads to a lower T e and less ionization. But for £ > 1000, a higher temperature 
luminosity produces a higher T e and more ionization. For example, at T = 10 8 K, the CNO 
elements are completely ionized at £ = 100, S becomes ionized at £ ~ 1000, and everything 
is ionized at £ ~ 10 4 . However, for T = 10 9 K, CNO do not become ionized until £ ~ 500, S 
still becomes ionized at £ = 1000, but everything is ionized at £ ~ 5000. These results can be 
understood since the ionization potential is oc Z 2 , where Z is the atomic number. The higher 
temperature emission has higher energy photons, which are less efficient at ionizing atoms 
with lower ionization potentials (compared to the photon energy), but are more effective at 
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ionizing atoms with high ionization potentials. 

For t w ~ 1, we expect that (3 ~ 0.5, and the value of £ mainly depends on v s . For u s4 > 1, 
the gas is completely ionized, while for v s $ = 0.5 only partial ionization is likely. Another 
issue is whether the preshock medium is optically thick in the photoionization continua, and 
the incident radiat i on fie ld is depleted in crossing the circumstellar gas. Using the results of 



Tarter fc Salpeterl (119691 ) . we find that at the higher temperatures of interest (~ 10 9 K) for 
the ionizing radiation, photon depletion is negligible and ionization remains nearly constant 
with radius. At lower temperatures for the radiation (~ 10 8 K), there is a transition to the 
Stromgren regime. 

The situation changes if the wind is moderately optically thick to electron scattering. 
An effect of electron scattering is to reverse the motion of photons, so that the energy density 
of ionizing photons is increased; the time spent by photons in a region is increased by ~ t w 
(equation IT] ), so that the rate of ionizations (and £) is increased by the same factor (e.g., 



Rossi 119791 ). While this effect favors higher photoionization and the escape of X-rays, there 
are several effects that disfavor high X-ray emission. First, Compton degradation in the wind 
leads to the loss of the high energy photons, as discussed above; emission above 2 keV is 
suppressed when t w ~ 16. Second, at moderate optical depths inverse Compton cooling tends 
to dominate bremsstrahlung (Fig. 1), so that X-ray emission is a smaller fraction of the shock 
power; inverse Compton is more important at higher shock velocities. In addition, initially 
outgoing X-ray photons have some chance of being scattered back and absorbed by the dense 
shell. Finally, at optical depths > 10, these two effects are likely to be larger than the increase 
in the ionizing radiation field, so that the ionization parameter is decreased and there is an 
increased chance of photoabsorption of the X-ray emission. Photoabsorption is important 
at low X-ray energies while Compton degradation is important at high X-ray energies. An 
accurate calculation of the X-ray emission is complicated and beyond the current paper. 
An ana logous physical situati on is X-ray emission from optically thick accretion onto white 



dwarfs fcvlafis &: Lamblll982l ). The photon trapping limit r w = c/v s here corresponds to an 



accretion rate at the Eddington limit in the accretion case. A difference is that the emission 
from the white dwarf surface can maintain complete ionization of the preshock gas, which is 
not necessarily the case here. 



3. IMPLICATIONS FOR OBSERVATIONS 

Most of the X-ray observations of Type Iln supernovae in dense media are an age > 1 
yr and sometimes much greater, so that the objects do not fall in the parameter space 



shown in Fig. 1. However, SN 2006gy, with an estimated « 10 ( jSmith McCrayl 120071 ; 
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Chevalier fc Irwinl 120111). was observed with Chandra on 2 006 N ov 14, at an age of 3 — 4 
months ( Smith et al.l 2007 ; Qfek et al. 2007 ). Smith et al. ( 2007 ) inferred a detection, with 



all counts at energies < 2 keV and an unabso rbed 0.5 — 2 keV lu minosity of 1.65 x 10 39 erg s 1 



assuming T = 1 keV. From the same data, lOfek et al.l (120071 ) inferred a nondetection with 



an upper limit of 1.6 x 10 erg s , assuming a photon index of 1.8. In either case, the X-ray 
l uminosity was mu ch less than the ob served photospheric lumin osity of ~ 3 x 10 44 erg s -1 
JSmith et al.l boioh . In the model of Ichevalier fc Irwinl J201lh for SN 2006gy, the shock 
wave radiation broke out in the mass loss region, so the optical depth outside the shock was 
initially c/v s . The para meters for th e X-ra y observation give -D*(i/10 day) -1 ~ 1, near the 



T w — c/v s line in Fig. 1. ISmith et al.l ( 120101 ) estimate v s ~ 4000 — 5000 km s 1 near the time 



of peak optical luminosity so, as discussed in Section 2.2, Comptonization by itself could 
limit the escaping photons to < 0.14 keV, and the loss of ionizing radiation would allow 
photoabsorption. 

Our discussion suggests the following sequence, as the density of surrounding mass loss 
increases. At low density, the optical luminosity is dominated by radioactivity and shock 
heating of the progenitor; X-ray emission from interaction is initially a small part of the 
luminosity, although it might become a larger part at late times when other power sources 
fade. As the optical depth in the wind approaches unity, interaction typically dominates the 
power input and X-ray emission can be a signficant part of the power. At higher densities, 
the X-ray emission falls relative to the optical luminosity because of inverse Compton cooling 
of the shocked region, Compton degradation in the wind, and photoabsorption. 
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Fig. 1. — Regimes of shock structure in terms of wind density/age vs. shock velocity. The 
lines in the figure mark the limits of a) presence of viscous shock (thick solid line), b) 
optical depth unity (dashed), c) rapid cooling (solid), d) cooling by inverse Compton vs. 
bremsstrahlung (solid), and e) electron-proton equilibration (dash-dot). 



